Interference in transport through double barriers in interacting quantum wires 
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We investigate interference effects of the backscattering current through a double-barrier structure 
in an interacting quantum wire attached to noninteracting leads. Depending on the interaction 
strength and the location of the barriers, the backscattering current exhibits different oscillation and 
scaling characteristics with the applied voltage in the strong and weak interaction cases. However, 
in both cases, the oscillation behaviors of the backscattering current are mainly determined by the 
quantum mechanical interference due to the existence of the double barriers. 
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I. INTRODUCTION 

As a fundamental many-body physical model, one- 
dimensional(lD) interacting electron systems are an ever- 
lasting research topic. Unlike its high-dimensional coun- 
terparts, which are well understood from the quasipar- 
ticle picture in the Fermi liquid theory,— ID interact- 
ing electron systems can be described by bosonic gapless 
collective excitations of the fermion density fluctuations 
within the framework of the Tomonaga-Luttinger liquid 
theoryj^"— 

Tomonaga-Luttinger liquid behaviors in ID systems 
have been revealed in measuring the transport proper- 
ties of cleaved-edge overgrowth quantum wires,— quan- 
tum Hall systems,^ and single-wall carbon nanotubesj^ 
It opens up the possibilities to test some theoretical pre- 
dictions even in simplified pure physical models, and trig- 
gers intensive attention to the physics of the Tomonaga- 
Luttinger liquids.* 

In mesoscopic Fermi-liquid systems, electronic trans- 
port is modeled by the transmission of an incident elec- 
tron through a potential barrier, and the conductance 
of the system is directly related to the transmission 
probability.^ However, potential barriers play a counter- 
intuitive role in the Tomonaga-Luttinger liquids. In a 
seminal paper, Kane and Fisher— have shown that the 
barrier is irrelevant for attractive electron-electron inter- 
actions, and cuts the system into two pieces for repulsive 
interactions. In transport measurements, Tomonaga- 
luttinger liquid quantum wires of finite length need to 
be connected to leads acting as electron reservoirs. If 
the reservoirs are modeled as ID noninteracting systems, 
then the model of inhomogeneous Tomonaga-luttinger 
liquids with different interaction strengths in different 
parts is appropriate to investigate the transport prop- 
erties of the Tomonaga-luttinger liquids. ^^'^^ It was ar- 
gued that the dc conductance is not renormalized by the 
interactions Jii^ An interesting phenomenon of Andreev- 
like reflections has been shown at the interfaces between 
the interacting wire and the noninteracting leads due to 
the mismatch of the interaction strengths. "'^^ The pres- 



ence of an impurity results in interference of the bosonic 
excitations which can be modulated by an applied bias 
voltage, and thus leads to characteristic oscillations of the 
backscattering current as a function of the dc voltage.— 
Feldman et alJ^ have showed that backscattering off a 
weak dynamic impurity would enhance the current. Pega 
et aliiS and Recher et alJ^ have demonstrated interesting 
Fabry-Perot interference patterns of the non-linear con- 
ductance as a function of the bias voltage in carbon nan- 
otubes attached to metallic reservoirs, where backscat- 
tering processes mainly occur at the two metal-contact- 
nanotube interfaces. 

In this paper, we investigate interference effects of elec- 
tron tunneling through double barriers in an interacting 
quantum wire. We assume that the ID interacting elec- 
tron system is adiabatically attached to two ID nonin- 
teracting electron systems, which act as electron reser- 
voirs. Such an assumption excludes the possibility of 
backscattering at the interfaces between the interacting 
and noninteracting segments. Backscattering events take 
place at the positions of the double barriers. Due to 
the existence of the Andreev-like reflections at the inter- 
faces, electron tunneling through the double barriers in a 
Tomonaga-Luttinger-liquid wire is more complicated and 
thus expected to be more interesting. As we show later, 
there exists a contribution from the quantum mechanical 
interference of bosonic excitations to the backscattering 
current, besides the contributions from backscattering off 
the separate barriers. It is the quantum mechanical in- 
terference that dominates the oscillatory characteristics 
of the backscattering current as a function of the ap- 
plied bias voltage. The aim of this work is to make clear 
how the oscillation pattern of the backscattering current 
in an interacting quantum wire with double barriers is 
modulated by the electron interaction strength and the 
arrangement of the double barriers. 



II. MODE AND FORMULATION 

We consider an interacting quantum wire of length L 
with double point-like barriers, which is at its ends con- 
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nected adiabatically to two semi-infinite noninteracting 
quantum wires acting as electron reservoirs with electro- 
chemical potentials iil and /i/j . The model Hamiltonian 
can be written as 



H = Hw + Hi 

f + OO 



2m dx^ 



^i(x)\'4}{x) 



+ U{x)^P{x)^?{x)Y 
Hb ~ I dxip\x)\\8{x — xi) + \5{x — X2)\i^{x), 



where Hw describes both the interacting wire and the 
interacting leads with varied electron-electron interac- 
tion strengths, i.e. U{x) = [/ if -L/2 < x < L/2 
and U{x) = elsewhere, Hb represents the Hamilto- 
nian for the point-like barriers, and the function ^{x) = 
jiL [1 — d{x -I- L/2)] -I- IjLb(){x — L/2) denotes the externally 
tunable electrochemical potential. 

Since what we are interested in is the low-temperature 
transport properties, it is convenient to reformulate the 
problem in the framework of standard bosonization.^^ 
First linearize the energy spectra about the Fermi points 
ikp, and introduce two species of electron operators 
ipa{x){a = R/L = ±) to describe the right-moving and 
left-moving fermions, then the excitations of the fermion 
system can be described by a bosonic field In this 

way the Hamiltonian for the fermion system can be recast 
in the following bosonic form 
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Acos(v47ri?(xi,i) -I- 2kFXi) 

+ Acos(V47r-!9(x2,t) + 2kFX2), 



g{x) lv{x) 

dxE{x,t)d{x,t), 



(1) 



(2) 



(3) 



where g{x) = [1 + U{x)/TThvF]~'^^'^ is the interaction pa- 
rameter, —eE{x,t) — dx^Ji — —^l5(x + L/2) + /i/?(5(a; — 
I//2), and v{x) = vp/gix) is the charge density wave 
velocity. 

In terms of the Bosonic field i9, the current operator 
can be written as 



i{x,t) 



=dt^ix,t). 



(4) 



The current-bias relations of the nonequilibrium sys- 
tem is of interest, it is appropriate to adopt the Keldysh 
formalism. Following the path integral technical proce- 
dures developed in Ref. 13 for a static impurity, the 
average current is obtained 



{I{x,t)) = {Io-lBs{x,t)). 



(5) 



In Eq. ([5]), (/q) = e^V/h is the background cur- 
rent in the absence of barriers, and (7^5(2;, t)) — 
_^^ d..^ijj^{x,t)4!R{x,t). -^ ^j^^ backscattering current. The 
average backscattering current is the sum of three con- 
tributions to leading order in A (see Appendix) 
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where wq = eV/h, Co{x,t; x' ,0) = (^■d{x,t)i}{x' ,0) 



[i?^(a;, t) + d'^{x' , 0)]/2y are the correlation functions of 

the bosonic field d{x,t) in the clean wire limit. These 
correlation functions can be obtained by expanding the 
Bosonic field on the basis of the eigenfunctions, which 
satisfy a specific inhomogeneous equationj^^ It is evi- 
dent that, besides the independent contributions from 
backscattering by the separate barriers, there exists a 
quantum mechanical interference contribution to the 
backscattering current due to the coexistence of double 
barriers. It is just such an interference term that gives 
rise to some interesting transport characteristics in inter- 
acting quantum wires with double barriers. 

At zero temperature, the correlation functions yield a 
simplified expression 
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where = Xi/L, a = w^/wc is the dimensionless cutoff, 
with LUi^ = Vp / gL being the inverse of the traversal time 
of the charge density wave and ojc being the high-energy 
cutoff frequency, ^ — {1 — g) / {1 + g) is the Andreev-like 
refiection coefficient at the interfaces. 

In terms of the dimensionless variables u = eV/hwL, 
T = ujLt, the backscattering current is rewritten into the 
following simple form in unit of e{\uj^j^ / uiuj^)"^ / fi^ uj l 
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The above expressions for the backscattering current 
Ibs ahow us to evaluate it for arbitrary values of the bar- 
rier positions (a;i, X2), the interaction strength g, and the 
applied voltage V . It is noticed that the backscattering 
current contributed from an individual barrier is depen- 
dent on the relative distance of the barrier to the wire 
center(|6/2l): while that from the quantum mechanical 
interference term is dependent on the location details of 
the double barriers(|6— and ^i-l-^). It seems impossi- 
ble to obtain an explicit expression for the backscattering 
current, we resort to numerical calculations to investigate 
the backscattering current in some typical cases. 



III. RESULTS AND DISCUSSIONS 
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First, we investigate the dependence of the backscat- 
tering current on the interaction strength. It is well 
known that the interaction strength is characterized by 
the interaction parameter 5 = [1 -f- U l'K'hvp\~^l'^ . We 
have < (7 < 1 for repulsive interactions and g — \ for 
noninteracting limit. It should be noted that the smaller 
the value of the stronger the interaction. In Fig. 1, 
we present numerical results of Ibs as a function of the 
applied voltage u for different typical interaction param- 
eters g ~ 0.25,0.5,0.75, as the barriers locate symmetri- 
cally near the ends of the wire with = ±0.4. Since 
the barriers are symmetrically located, the backscatter- 
ing currents arising from different barriers are the same. 
It can be observed that, as the interaction decreases, 
the oscillation of the incoherent addition of the backscat- 
tering currents contributed from the two barriers disap- 
pears gradually, while that contributed from the quan- 
tum mechanical interference persists and is more promi- 
nent. The oscillatory behavior of the backscattering cur- 
rent is mainly determined by the quantum mechanical 
interference of the bosonic waves backscattering off the 
double barriers. We attribute the less pronounced oscil- 
lation of the backscattering current in the strong inter- 
action case to the suppression of the quantum mechani- 
cal interference by the electron-electron interaction. As 
shown by Dolcini et al.^'^, the oscillation of the backscat- 
tering current in the single barrier case arises from a 
combined effect of the barrier, the finite length and the 
interaction in the wire. The phase shift of bosonic exci- 
tations traveling between the barrier and the interfaces 
is responsible for such an oscillation and can be modu- 
lated by an applied bias voltage. In our case with double 
barriers, competition between the single-barrier interfer- 
ence and the double-barrier quantum mechanical inter- 



FIG. 1: Backscattering current Ibs (in unit of 
e{\Lu^j^/u)^)'^ /U^lul) as a function oi u = eV/hujL with 
the double barrier location = X1/2/L — ±0.4 for 

different interaction parameters (a) g = 0.25, (b) g — 0.5 and 
(c) g — 0.75, respectively. The solid line refers to the total 
backscattering current Ibs, the dot line to the incoherent 
addition of independent contributions from different barriers 
Ibs + ^BS ^-nd the dash line to the quantum mechanical 
interference term Ij^s- 



ference results in an interesting oscillatory behavior of 
the backscattering current. In the cases of a given bar- 
rier location, the oscillation period of the backscattering 
current contributed from the quantum mechanical inter- 
ference is irrespective of the interaction strength, while 
the oscillation period from the incoherent addition of the 
backscattering currents off different barriers is increased 
as the interaction parameter is decreased, and eventually 
becomes infinite in the noninteracting limit, which can 
be found in Fig. 1. 

Then we consider two special cases for strong interac- 
tion {g — 0.25) where the double barriers locate symmet- 
rically near the midpoint of the interacting wire (6/2 = 
±0.01) and right at the ends of the wire {£,1/2 = ±0.5). It 
can be found from Fig. 2 that in both cases the backscat- 
tering current oscillates in a more pronounced way, com- 
pared to the case where the double barriers are located 
symmetrically near the interfaces {£,1/2 — ±0.4, Fig. 1 
(a)). In Fig. 2, we also observe that, the oscillation fre- 
quency of the incoherent addition of the backscattering 
currents off the separate barriers is twice that due to the 
quantum mechanical interference when the double bar- 
riers are located at the ends of the wire, while the two 
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FIG. 2: Backscattering current Ibs (in unit of 
e^Xui^^/ujcY /Ti^ujl) as a function of u = eV/TiuJL as 
g = 0.25 with different double barrier locations (a) 
(.1/2 = = ±0.01 and (b) ^1/2 = = ±0.5. 

The solid line refers to the total backscattering current 
Ibs, the dot line to the incoherent addition of independent 
contributions from different barriers I^g + Ibs ^-^d the dash 
line to the quantum mechanical interference term Ij^s- 



frequencies are about the same when the double barriers 
are located near the midpoint of the wire. The reason 
is that the round-trip ballistic time for bosonic excita- 
tions propagating between a barrier and an interface in 
the former case is twice that in the later case. 

Up to the present we have analyzed the current-voltage 
characteristics in cases where the two barriers are sym- 
metrically located. As the barriers are located asymmet- 
rically, it is expected that such an asymmetry has trivial 
effects on the current- voltage characteristics of weakly in- 
teracting wires, but influences significantly the transport 
properties of strongly interacting wires. In Fig. 3, we 
provide the results of the backscattering current of the 
interacting wire with symmetric and asymmetric loca- 
tions of the barriers in both weak interaction((7 = 0.75) 
and strong interaction(g = 0.25) cases. We find that 
in the weak interaction case, the backscattering currents 
are about the same for symmetric and asymmetric bar- 
rier locations(Figs. 3(a) and 3(b)), as long as the spacing 
between the double barriers is fixed. However, such a sce- 
nario is changed in the strong interaction case, where the 
current oscillation strongly depend on the symmetry of 
the barrier locations (Figs. 3(c) and 3(d)). It is interest- 
ing to note that, the period of current oscillation depends 
strongly on the spacing between the double barriers in the 
weak interaction limit, i.e., the bigger the barrier spac- 
ing, the smaller the oscillation period. This phenomenon 
can not be observed in the strong interaction limit. We 
attribute the dependence of the current oscillation on 
the barrier spacing in the weak interaction limit to the 
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FIG. 3: Backscattering current Ibs (in unit of 
e{\ujf^ / Lj^)'^ /H^lul) as a function oi u = eV/hcuL with 
different double barrier locations and for different inter- 
action parameters (a) g — 0.75, ^1/2 ~ X1/2/L — ±0.2 
(solid), Ci/2 = X1/2/L = 0,0.4 (dotted) (b) g = 0.75, 
Ci/2 = a;i/2/L = ±0.3 (solid), C1/2 = S1/2/L = -0.15,0.45 
(dotted), (c) g = 0.25, ^1/2 = xi/2/L = ±0.2 (solid), 
(1/2 = x-i_/2/L = 0,0.4 (dotted), and (d) g = 0.25, 
(1/2 = xy2/L = ±0.3 (solid), Ci/2 = S1/2/L = -0.15,0.45 
(dotted). 



resonant tunneling through a double-barrier structure in 
one-dimensional noninteracting electronic systems*'^^ In 
the case of strong interactions, different locations of the 
barriers lead to different interference patterns for bosonic 
plamonic excitations traveling between a barrier and the 
interfaces, and then results in different oscillation behav- 
iors of the backscattering currents contributed from the 
separate barriers. 

Finally, we would like to point out that, depending on 
the interaction strength, the backscattering current ex- 
hibits different scaling rules with the applied bias voltage 
in a finite-length Tomonaga-Luttinger-liquid wire with 
double barriers, as can be found from Figs. I, 2, 3. This 
observation deserves further investigation. 



IV. CONCLUSION 

In summary, we have investigated interference effects 
in electron transport through a double-barrier structure 
in an interacting quantum wire, which is ideally attached 
to two noninteracting leads. We have found that the os- 
cillation of the backscattering current with the applied 
voltage is mainly determined by the quantum mechanical 
interference due to the coexistence of the double barri- 
ers. It is contrast to the single barrier case, where the 
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current oscillation is strongly dependent on the interfer- 
ence of the Andreev-like reflected plasmonic excitations 
propagating between the single barrier and the interfaces. 
Depending on the interaction strength and the locations 
of the barriers, the competition between these two kinds 
of interferences results in different oscillation and scal- 
ing behaviors of the backscattering current with the bias 
voltage. 
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Appendix A: Outline of derivation of backscattering 
current expression 

In this appendix,we outline the main formulae to cal- 
culate the backscattering current based on the Keldysh 
functional approach following Ref. 13. 

The generating functional takes the form 



^S + ^l dxJ(x)^(x) ,(A1) 



where the action functional of the system can be written 
in terms of the boson field 'd{x,t) as 



h dx dt 



1 



1 



2g{x) lv{x) 



dx / dtE{x,t)-d{x,t) - / dtHB- (A2) 



V2 



After introduce the standard keldysh time contour and 
denote by i?+ and the complex fields on the upper 
and lower time branches of the Keldysh contour, and 
define four Green's functions averaged with respect to 

the free Hamiltonian Cj^''' {r' ; r" ) = (z9''(r')t?'''(r"))o, the 
generating functional is rewritten as 



Z[J] = ^ j Dd^e^v{-\j dv'dv" z9"(r') 

j; j;'=± 

(C-i)W(r';r")^"'(r")}exp{ ^ j dv'E{v') 



+00 



dt'EBl^"^] 



7,=± 

i 

V2 



dxJ(x)^''(x)))}, 



where (r; r') is the inverse of of a 2 x 2 matrix formed 
from the above four Green's functions. 
Define the following matrices 



i9 




and 



Co 



C7o++(r,r') Co+'(r,/) 



(A3) 



(A4) 



the generating functional is reduced to the following sim- 
plified form 



Z[J] 



— / ,D^e(-^"C-^i9+2a"Qi9)/2 



X exp { - J ^ 77 / dt'i7B[i?'']}.(A5) 

■q=± 

After Shift the fields i? ^ + Aj, Aj = iCoQ'^J, one 
obtains a factorized form of the generating functional 

Z[J] = e-J"QCoQ"J/2 

'exp(-j5]77 / HB[r + A'}]dt')) .(A6) 



From the expression of current operator Eq. we 
have 

{I{x,t)) = ^dt{d{x,t)) 



= 1 dr'Co^(x;r')i?(r')--^ 



—d \— 

dt 



V2h 



(E 

e^V/h - Ibs 



(5J(x) 



(A7) 



where C^(r;r') = e{t ~ t')([7?(r), z9(r')])o is the retarded 
Green's function, and the backscattering current takes 
the following form 



IBS 



^ ^ /dt'ao(x;r00-+(r',))^.(A8) 



2=1,2 



Here ao(r;r') = 2ie^9tC(f (r; r')//i is the local con- 
ductivity of clean wire, and the backscattering cur- 
rent operator j^(x) = -ji^^^^EW^T^ ^^i^^ ^o('') = 
J dx'C(f (r; x')£^(x'), (•••)_> denotes an average along 
the Keldysh contour with respect to the shifted Hamil- 
tonian = Ho[S] +Hb[^ + Ao]. 

Some algebras after substitution of the expression of 
Hb into Eq. (|X8)) finally yield 
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ain[V^'&^ixi,t') + 2kFX,+ujQt']) . (A9) 
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